Super-resolution single-beam imaging via compressive sampling 
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Based on compressive sampling techniques and short exposure imaging, super-resolution imaging with thermal 
light is experimentally demonstrated exploiting the sparse prior property of images for standard conventional 
imaging system. Differences between super-resolution imaging demonstrated in this letter and super-resolution 
ghost imaging via compressive sampling (arXiv. Quant-ph/0911.4750vl (2009)), and methods to further 
improve the imaging quality are also discussed. © 2010 Optical Society of America 

OCIS codes: (100.6640) Superresolution; (110.6150) Speckle imaging; (100.3010) Image reconstruction techniques. 
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Back in the 19th century, Lord Rayleigh had demon- 
strated that the best resolution ol conventional imag- 
ing system was determined by the numerical aperture 
of the imaging lens [1]. Some decades later, the true 
limit on imaging arose from the optical wavelength A 
and the recoverable resolution was A/2 exploiting the 
evanescent components of the electromagnetic held near 
the surface of object [2-5] . Several digital image process- 
ing approaches can realize super-resolution, but they are 
highly sensitive to noise in the measured data and to 
the assumptions made on the prior knowledge [6-9]. In 
recent ten years or so, a new imaging technique called 
ghost imaging (quantum imaging) can realize subwave- 
length interference with the resolution beyond the clas- 
sical diffraction limit by a factor of 2 using entangled 
source or thermal light source [10-12]. When signals sat- 
isfied certain sparsity constraints, Donoho has demon- 
strated mathematically that super-resolution restora- 
tion is possible below the Rayleigh threshold (Q -C 3, 
where A is the lattice span and $7 is the frequency cut- 
off) [13] and a new information processing techniques 
called as compressive sampling (CS) provides a solver for 
large-scale sparse reconstruction [14-18]. Combining the 
sparse prior property of images with the ghost imaging 
method, super-resolution ghost imaging via compressive 
sampling (GICS) with the resolution beyond the diffrac- 
tion limit by a factor of 10 has been demonstrated ex- 
perimentally even without looking at the object, and 
its physical principle suggests that using compressive 
sampling reconstruction algorithms, conventional imag- 
ing system may also realize super-resolution because all 
objects have sparse representations when expressed in 
the proper representation basis [19]. For the conventional 
imaging system widely applied in people's life, as shown 
in Fig. 1, can we directly obtain super- resolution imaging 
with thermal light exploiting the sparse prior property 
of images and CS reconstruction algorithms? 

Here, we demonstrate super-resolution imaging via 




glass Object 



Fig. 1. Standard schematic of conventional imaging and 
super-resolution imaging via compressive sampling with 
thermal light. 



compressive sampling borne on thermal light using the 
standard conventional imaging system shown in Fig. 1. 
The pseudo-thermal light source S, which is obtained by 
passing a laser beam (with the wavelength A=532nm and 
the source's transverse size fl=4mm) through a slowly 
rotating ground glass disk, hrst goes through an object 
and then to an imaging lens with the focal length /. 
The intensity distributions transmitted through the ob- 
ject are finally recorded by a CCD camera. The distances 
z\ , z-i and the focal length of the lens / obey the Gaus- 
sian thin- lens equation: — + — = i According to the 
Rayleigh criterion [1] , the resolution limit Ax of conven- 
tional imaging is determined by the wavelength A and 
the numerical aperture (N.A.) of the lens /, namely 



Ax = 0.61 
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where L is the effective transmission aperture of the 
imaging lens /. 

Based on the theory of statistical optics [20], the in- 
tensity distribution on the detection plane at a certain 
time s can be expressed as by the Fresncl diffraction in- 
tegral [7] 



I s (x,y)= I dxodyodx' dy' E s (x ,yo)E*(x' ,y' Q ) 

xh*(x,y;x' ,y' )h(x,y;x ,y ). (2) 



where E s (xo,yo) denotes the light field on the source 
plane at time s, h(x, y; Xq, yo) and h*(x, y; x' , y' ) are the 
impulse function of optical system and its phase conju- 
gate, respectively. 

For the specific optical system depicted in Fig. 1 and 
under the paraxial approximation, the impulse response 
function of the optical system is 

h(x,y;x ,yo) oc / dx'dy' exp{ — (x 2 + y' 2 )} 
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where T(x, y) is the transmission function of the object, 
sin c(x) = s and sin c(y) = ain ^ y > . 
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Similar to the idea of ghost imaging via compressive 
sampling (GICS) [19], for the schematic shown in Fig. 
1, the CCD camera has recorded the similar but not the 
same information among each observation because of the 
fluctuation of the light field, thus the intensity distribu- 
tions recorded with the CCD camera obey Gaussian sta- 
tistical distribution and also satisfy the requirement of 
incoherence measurements of CS. Different from GICS, 
both the sensing basis and the information transmitted 
through the object are recorded by the same CCD cam- 
era. If the intensity distributions located at the CCD 
camera plane at time s are described by I s (x,y), then 
the image of the object can be reconstructed by solving 
the following convex optimization program [17]: 

\T CS \ = \T'\; subject to: 
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(4) 



where K is the acquisition numbers, r is a nonnega- 
tive parameter, ||T^|| 2 denotes the Euclidean norm of V 
||y|| 1 = J2 i \vi\ is the l\ norm of V, and \T C s\ is the ob- 
ject's function recovered by CS reconstruction algorithm 
Moreover, I s (x,y) in Eq. (4) is 
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Fig. 2. Experimental results of conventional imaging and 
super-resolution imaging via compressive sampling of 
two objects using the schematic shown in Fig. 1. Left 
column: a double-slit (slit width a=90//m and center-to- 
ccnter separation d=180/imm); right column: an aper- 
ture ("zhong" ring), (a). The original object; (b) and 
(c) are the images reconstructed by conventional imag- 
ing and super-resolution imaging via compressive sam- 
pling with thermal light, respectively (with 1000 obser- 
vations). 



In the experiment, the distances listed in Fig. 1 were as 
follows: z=200mm, zi=Z2=800mm, and the focal length 
of the lens /=400mm. The effective transmission aper- 
ture of the imaging lens was L=2mm and the exposure 
time widow for the CCD camera was set to be 1 ms. 
The experimental results of imaging a double-slit and 
an aperture ("zhong" ring) are shown in Fig. 2. For 
comparison, Fig. 2(a) displays the images of the objects 
by an ordinary lens imaging system using thermal light. 
The images of the objects for conventional imaging, as 
depicted in Fig. 2(b), appear in the screen of the CCD 
camera. Using the method of imaging via compressive 
sampling and the gradient projection for sparse recon- 
struction (GPSR) algorithm [17], we observe the super- 
resolved images of the objects in Fig. 2(c). 

The differences of super-resolution imaging demon- 
strated here and super-resolution ghost imaging shown 
in Ref. [19] greatly deserve discussion. GICS is two-arm 
nonlocally imaging system and the diffraction limit of 
optical system is determined by the transverse coher- 
ence length on the object plane in the test path. In 
the image extraction process of CS reconstruction al- 



gorithm, the sensing basis in CS, highly correlated with 
the light field on the object plane in the test path, is reg- 
istered by the CCD camera in the reference path and a 
bucket detector fixed in the test path has received all the 
information transmitted through the object. However, 
for the standard conventional imaging system shown in 
Fig. 1, because of the diffraction effect caused by the fi- 
nite transverse size of the lens, the CCD camera records 
spatial low-resolution intensity distributions of the light 
field transmitting through the object and they are used 
to compose the sensing basis of CS reconstruction tech- 
niques. Apparently, the sensing basis of this two super- 
resolution imaging methods is different. For the former, 
the intensity distributions outside the field of the ob- 
ject also are incoherent among each observation and are 
helpful to the reconstruction of images. While for the 
latter, only the intensity distributions within the field 
of the object are useful to the measurements because 
the intensity distributions outside the field of the object 
are always zero for each observation (namely completely 
coherent). Therefore, the vector length of the represen- 
tation basis for GICS will be much larger than that of 
conventional imaging system because it is the same as 
the vector length of sensing basis in CS reconstruction 
algorithms. Based on the characteristic of CS reconstruc- 
tion theory, the larger the sparse degree of images in 
the representation basis is, the higher the probability 
of stable recovering the object will be [14-16]. Hence, 
GICS will have much higher ability to obtain supcrrc- 
solvcd images than that in conventional imaging system 
because the images of the same objects are much sparser 
in the representation basis with larger vector length. In 
addition, from the viewpoint of speckle, for the above 
two super-solution imaging methods with thermal light, 
super-resolution imaging relies on multiple short expo- 
sures and makes full use of the similar but not the same 
information among each observation. Based on the basic 
principle of speckle imaging [20] , optical imaging meas- 
ured over an exposure time on the order of several sec- 
onds (namely long exposure imaging) will give an aver- 
age image and can not observe the intensity fluctuation 
of the light field. Completely different from long exposure 
imaging, short exposure imaging gives a speckled appear- 
ance and of course can catch the fluctuation of the light 
field. And the shorter the detection time is, the stronger 
the intensity fluctuation of the light field is and help- 
ful to recovering the objects' images. Further, combining 
the improved images obtained by previous digital image 
processing approaches (such as dcconvolution, analytical 
continuation) with suitable compressive sampling recon- 
struction algorithms, the ability of super-resolution may 
be further enhanced. 

In conclusion, combining the sparse representation of 
the images in the proper representation basis with suit- 
able compressive sampling reconstruction algorithms, su- 
perresolved images can also be obtained in standard con- 
ventional optical imaging system. This imaging skill can 
be directly applied to people's life and is very useful to 



the long-range imaging and microscopy. 
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